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Preheating in the Standard Model with the Higgs-Inflaton coupled to gravity 
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We study the details of preheating in an inflationary scenario in which the Standard Model Higgs, 
' strongly non-minimally coupled to gravity, plays the role of the inflaton. We find that the Universe 

, does not reheat immediately through perturbative decays, but rather initiate a complex process in 

■ which perturbative and non-perturbative effects are mixed. The Higgs condensate starts oscillating 

' around the minimum of its potential, producing W and Z gauge bosons non-perturbatively, due 

to violation of the so called adiabaticity condition. However, during each semi-oscillation, the 
created gauge bosons partially decay (perturbatively) into fermions. The decay of the gauge bosons 
prevents the development of parametric resonance, since bosons cannot accumulate significantly at 
the beginning. However, the energy transferred to the decay products of the bosons is not enough 
to reheat the Universe, so after about a hundred oscillations, the resonance effects will eventually 
dominate over the perturbative decays. Around the same time (or slightly earlier), backreaction 
from the gauge bosons into the Higgs condensate will also start to be significant. Soon afterwards, 
(~| , the Universe is filled with the remnant condensate of the Higgs and a non-thermal distribution 

^ of fermions and bosons (those of the SM), which redshift as radiation and matter, respectively. 

I ' We compute the distribution of the energy budget among all the species present at the time of 

backreaction. From there on until thermalization, the evolution of the system is highly non-linear 
p:^ ■ and non-perturbative, and will require a careful study via numerical simulations. 
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I. INTRODUCTION 
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\^ . Inflation is nowadays a well established paradigm, consistent with all the observations, that solves most of the 
Tj" ■ puzzles of the Hot Big Bang Model in a very simple and elegant way. It is able to explain not only the homogeneity 
, and isotropy of the present Universe on large scales, but also the generation of almost scale invariant primordial 
y—{ perturbations that give rise to the structure formation However, the naturalness of inflation is directly related 
00 i to the origin of the inflaton. Most of the inflationary models proposed so far require the introduction of new degrees 
■ of freedom to drive inflation. The nature of the inflaton is completely unknown, and its role could be played by any 
^ ' candidate able to imitate a scalar condensate (tipically in the slow-roll regime), such as a fundamental scalar field, 
a fermionic or vector condensate, or even higher order terms of the curvature invariants. The number of particle 
physics motivated candidates is as big as the number of extensions of the Standard Model (Grand Unified Theories, 
supersymmetry, extra dimensions, etc.), where it is not very difficult to find a field that could play the role of the 
infiaton 

In addition, given a model we must find a graceful exit to infiation and a mechanism to bring the Universe from a 
cold and empty post-inflationary state to the highly entropic and thermal Friedmann Universe [3| . Unfortunately, the 
theory of reheating is also far from being complete, since not only the details, but even the overall picture, depend 
crucially on the different microphysical models. It seems difficult to study the details of reheating in each concrete 
model without the experimental knowledge of the strength of the interactions among the infiaton and the matter 
fields. Because of this, most of the work until now has focussed on models encoding the different mechanisms that 
could play a role in the process, with the strength of the couplings set essentially by hand. The relative importance of 
each one of these mechanisms can only be clarified in light of an underlying particle physics model, able to provide us 
with the couplings among the infiaton and matter fields. From this point of view it is very difficult to single out a given 
model of inflation, and even more difficult to understand the details of the reheating process via the experimental 
access to the couplings. 



* E-mail: juan.garciabellido@uam.es 
t E-mail: daniel.figueroa@uam.es 
t E-mail: javier.rubio@uam.es 



2 



We may be very far away from understanding the microphysical mechanisms responsible for inflation, but maybe 
the natural candidate for being the inflaton was already there long ago. If we do not want to introduce new dynamical 
degrees of freedom in the theory, apart from those present in the Standard Model (SM) of particle physics, and at the 
same time we require Lorentz and gauge invariance, we are left with just one possibility: the Higgs field. Early models 
of inflation in terms of a Higgs-like scalar fleld h with a quartic self-interaction potential jh"^ need an extremely small 
coupling constant A ~ 10~^^ f^, and are also nowadays excluded at around 3a by the present observational data [5]. 
However, the SM-based inflation may be rescued from these difficulties replacing the usual Einstcin-Hilbert action by 
a non- minimal coupling of the Higgs field to the Ricci scalar [1, [1, HjBH' ^ scalar-tensor theory fashion 

SiG = J d^xV^^H^HR. (1) 

The induced gravity (IG) action is indeed the most natural generalization of the Standard Model in a curved spacetime, 
given the relevance of the non- minimal coupling to gravity for renormalizing the theor y llOl . It belongs to a group 
of theories known as scalar-tensor theories, originally introduced by Brans and Dicke [11| to explain the origin of 
the masses. In those theories not only the active masses but also the gravitational constant G , are determined by 
the distribution of matter and energy throughout the Universe, in a clear connection with the Mach principle. The 
interaction that gives rise to the masses should be the gravitational one, since gravity couples to all particles, i.e. to 
their masses or energies. The gravitational constant G in these theories is replaced by a scalar function that permeates 
the whole space-time and interacts with all the ordinary matter content, determining how the later moves through 
space and time. Any measurements of an object's mass depends therefore on the local value of this new field. 

The successful Higgs mechanism lies precisely in the same direction of the original Mach's idea of producing mass 
by a gravitational-like interaction. The role of the Higgs field in the Standard Model is basically to provide the inertial 
mass of all matter fields through the local Spontaneous Symmetry Breaking mechanism^. The Higgs boson couples 
to all the particles in the Standard model in a very specific way , with a strengh proportional to their masses p3] , and 
mediates a scalar gravitational interaction of Yukawa type [IJ, [lj| , between those particles which become massive as 
a consequence of the local Spontaneous Symmetry Breaking. According to the Equivalence Principle, it seems natural 
to identify the gravitational and particle physics approaches to the origin of the masses. From this point of view, the 
induced gravity action ([1]) would be an indication of a connection between the Higgs, gravity and inertia. Indeed, the 
action ^ is, at least at the classical level, just a different representation of the Starobinsky's model of inflation p].[l5|. 
where inflation is entirely a property of the gravitational sector. Both representations of the same theory are simply 
related by a Legendre transformation. This fact, together with the possibility of having an inflationary expansion of 
the Universe, makes the model extremely appealing. Unfortunately, the induced gravity model cannot be accepted as 
a completely satisfactory inflationary scenario, since the gauge bosons acquire a constant mass in the Einstein frame 
and totally decouple from the Higgs-inflaton field which translates into an inefficient reheating of the Universe. 
Notice however that the action ([1]) is not the most general one that can be written in a nontrivial background. As 
was shown in [l6j the simultaneous existence of a reduced bare Planck mass Alp and a non-minimal coupling of a 
Symmetry Breaking field to the scalar curvature 

Shg = I d^x^l + iH^HR \ , (2) 
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avoid the decoupling of the gauge bosons and can give rise to an inflationary expansion of the Universe together with 
a potentially successful reheating. 

In this paper we initiate the study of the reheating process in the model presented in Ref. [16], in which the 
Symmetry Breaking field is the Standard Model Higgs, strongly non-minimally coupled to gravity and playing the 
role of the inflaton. The novelty and great advantage of this model is its connection with a well-known microphysical 
mechanism, hopefully accessible in the near future accelerator experiments. The measurement of the Higgs mass 
will complete the list of the couplings of the Standard Model and, therefore, one should be able to study all the 
details of the reheating mechanism. This makes the model under consideration extremely interesting and, potentially, 
predictive. Reheating in the context of scalar-tensor theories has been studied, in the Hartree approximation by 
Ref. ^17], and perturbatively by Ref. fisj, but without the Higgs boson playing the role of the inflaton, and therefore 
without an explicit coupling of the fundamental scalar field to matter. 

Studying the details of reheating in this Higgs-Inflaton scenario, we have put special attention to the relative impact 
of the different mechanisms that can take place. We have found that the Universe does not reheat immediately through 
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perturbative decays, but rather initiate a complex procces in which perturbative and non-perturbative effects are 
mixed. The Higgs condensate starts oscihating around the minimum of its potential, producing Z and W gauge bosons 
due to violation of the so called adiabaticity condition. During each semi-oscillation, the non-perturbatively created 
gauge bosons decay (perturbatively) into fermions. This decay prevents the development of the usual parametric 
resonance, since bosons do not accumulate significantly at the beginning. The energy transferred to the decay 
products of the bosons is not enough to reheat the Universe within a few oscillations, and therefore the resonance 
effects will eventually dominate over the perturbative decays. Around the same time, the backreaction from the gauge 
bosons into the Higgs condensate will also start to be significant. Soon afterwards, the Universe is filled with the 
remnant condensate of the Higgs and a non-thermal distribution of fermions and bosons (those of the SM), which 
redshift as radiation and matter, respectively. We end the paper computing the distribution of the energy budget 
among all the species present at the time of backreaction. From there on until thermalization, the evolution of the 
system is highly non-linear and non-perturbative, and will require a careful study via numerical simulations, to be 
described in a future publication. 

The paper is organized as follows. In section [IT] we present the model with the Higgs field non-minimally coupled 
to the scalar curvature, transform it into a new frame where the action takes the usual Einstein-Hilbert form, and 
derive an approximate inflationary potential. In section [TTIl we study the effect of the conformal transformation in the 
matter sector, including the interaction among the Higgs, vector bosons and fermions. Section llVI is devoted to the 
analysis of the different reheating mechanisms, both perturbative and non-pcrturbativc, that can take place, leaving 
for Section |V] the analysis of the combined effect of parametric resonance and perturbative decays. We then study 
the backreaction of the produced particles on the Higgs oscillations and the end of preheating in Section IVII Finally 
the conclusions are presented in Section IVlTl 



II. THE STANDARD MODEL HIGGS AS THE INFLATON 



The Glashow-Weinberg-Salam [12| action is divided into four parts: a fermion sector (F) which includes the kinetic 
terms for the fermions and their interaction with the gauge bosons, a gauge sector (G), including the kinetic terms 
for the intermediate bosons as well as the gauge fixing and Faddeev-Popov terms, a Spontaneous Symmetry Breaking 
sector (SSB), with a Higgs potential and the kinetic term for the Higgs field including its interaction with the gauge 
fields, and finally, a Yukawa sector (Y), with the interaction among the Higgs and the fermions of the Standard Model, 

SsM — Sp + Sg + SssB + Sy ■ (3) 

The simplest versions of this Lagrangian in curved spacetime follow the principles of general covariance and locality for 
both matter and gravitational sectors. To preserve the fundamental features of the original theory in flat space-time, 
one must also require the gauge invariance and other symmetries in flat space-time to hold for the curved space-time 
theory. The number of possible terms in the action is unbounded even in this case and some additional restrictions 
are needed. A natural requirement could be renormalizability and simplicity. Following this three principles (locality, 
covariance and restricted dimension), and the previously motivated requirement of not introducing new dynamical 
degrees of freedom, the form of the action is fixed, except for the values of some new parameters to be determined by 
the physics. This procedure leads to the non-minimal Lagrangian for the Standard Model in the presence of gravity, 
given by 

•S'sMG = SsM + Shg 1 (4) 

where ^sm is the Standard Model part ([3]) defined above, and Srg is the new Higgs-gravity sector, given by Eq. 
([2]). Here Mp = (SttG)"^/^ is the reduced Planck mass, R the Ricci scalar, H the Higgs field, and ^ is the announced 
non-minimal coupling constant. As showed in Ref. ^16;], the parameters ^ and the self-coupling A of the Higgs potential 
are related by ^ ~ 49000-\/A. In the unitary gauge, H = and neglecting all gauge interactions for the time 

being, the Lagrangian for the Higgs-gravity sector in the so-called Jordan (J) frame takes the form 



Shg + SssB 3 
where f(h) = {Ml + ih'^)/2, and 



.f{h)R--9^''d^hd,h^U{h) , (5) 



U{h) = \{h'-v^)\ (6) 



is the usual Higgs potential of the Standard Model, with vev v = 246 GeV. 
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In order to get rid of the non-minimal coupling to gravity, we proceed as usual, performing a conformal transfor- 
mation UM 



such that we obtain the Lagrangian in the so-called Einstein (E) frame 



^HG + ^SBS ^ / d,^^^ 



R + 3g'"'\7f,S/^lnn'^ - -g'^'^V^ In In f^^ 



(7) 



^^-^C/W • (8) 



The usual Einstein-Hilbert term can then be obtained imposing f{h)/il'^ = Mp/2, which implies the following relation 
between the conformal transformation and the Higgs field 



(9) 



This allows us to write the Lagrangian ([5]) completely in terms of h 



'HG "T '-'SSB 



M' 



f74 



g^-d^hd,h--^Uih) 



(10) 



where we have neglected a total derivative that does not contribute to the equations of motion. As we will be working 
in the Einstein frame from now on, we will skip over the tilde in all the variables to simplify the notation. 

Notice that the conformal transformation ([7|) leads to a non-minimal kinetic term for the Higgs field, which can be 
reduced to a canonical one by making the transformation 



dh 



dx _ n'^ + 6^^h'^/M^ _ l+^{l + 6()h^/Mj, 



(11) 



where x is a new scalar field. Doing this, the total action in the Einstein frame, without taking into account the gauge 
interactions, is simply 
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HG ^SSB 



/\f2 1 



(12) 



with 



v{x) 



u{Kx)) , 



(13) 



the potential in terms of the new field x- To find the explicit form of the potential in this new variable x, we must 
find the expresion of h in terms of x- This can be done by integrating Eq. pip , whose general solution is given by 



^X{h) = v/l + eCsinh-i ( Vl + 6e«) - v^sinh 



(14) 



where u = y/^h/Mp. Since ^ 1, we can take 1 -I- 6^ « 6^ and, using the identity sinh ^ x — ln{x + \^ x^ + 1) for 
— oo < x < CX3, we can approximate Eq. (I14p by 



^x(/^)-y^ln(l + .Y/^ 



or, equivalently. 



(15) 



(16) 



where a — ■y/2/3 and k = Mp^. The x field is therefore directly related in this aproximation (just in the limit ^ ^ 1 
and far from u = 0) to the conformal transformation SI in a very simple way and the inflationary potential (|13[) is 
just given by 



■Ml 



(17) 
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Mp 

FIG. 1: Comparative plot of the exact solution (red continuous line) obtained parametrically from Ea. (|14|) . the analytic formula 
P8[) for the potential (blue dashed line), and their parametrization (|20p (green dotted line). 



Since v ^ Afp, then 1 + i^jr ~ 1 and we can savely ignore the vev for the evolution during inflation and preheating, 
and simply consider the potential 

nx) = ^(l-e— (18) 

Notice that the previous potential only parametrizes partially the original potential ([S]), since it neglects the region 
X < 0, as can be seen in Figll] where we compare the exact solution (red continuous line) obtained parametrically 
from Eq. with the analitic formula (blue dashed line). Both solutions agree very well in the region of 

positive X but differ substancially for x < 0- The conformal transformation is even ill-defined in the negative field 
region. From Eq. ^ and Eq. (|16p we have 

^^^■^-1 = e^'^x „! = (!_ ^-'^'^xyo.'^x ^ (19) 

which is inconsistent, since the left-hand side of this equation is positive definite, while the right hand is negative 
definite for x < 0. Taking this into account, in order to study the different mechanisms of the post-inflationary regime, 
we will then use the parametrization 

yix) = ^{l-e—\^\)" . (20) 

which correctly describes the potential obtained from Eq. (HH), for the whole field range of interest. In Fig. [U this 
parametrization (green dotted line) is again compared to the exact solution (red continuous line) obtained from Eq. 
(|14p . Around then minimum, the potential (j20[) can be approximated as 

V{x)^\m\^ + ^V{x), (21) 

where = XMp/3£,^ is the typical frequency of oscillation and AV are some corrections to the quadratic ap- 
proximation, which soon become negligible after inflation ends, see section HVl Note, nevertheless, that approxima- 
tions (jlSp . p^ and therefore parametrization (PO)) . do not describe correctly the potential for very small values of the 
field u ^ X ^ Mp/^. As can be read from Eq. pT|) . for |x| ^ Xt = Mp/^, we have ^ ~ 1 and therefore, there is a 
transition in the potential from ((2T|) to V{x) ~ jX^- However, as will be shown in section [IV CI the transition region, 
Ixl < Xti is several orders of magnitude smaller than the non-adiabaticity region, |x| < Xa, see Eq. (I79P . inside which 
the concept of particle is not properly defined. Therefore, from now on we will neglect the change in the behaviour of 
the potential (from ^M^x^ to jx"^) i^i this "small" field region, since Xt ^ Xa- See section HVl for more details. 
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The analysis of the inflationary potential pp)) can be performed either in the Jordan [2^ [2l| or in the Einstein 
frame [l^ with the same result. The slow roll parameters can be expressed analytically, in the limit of 3> Mp/^ ^ 
t!^, as a function of x, 

_Ml(V'{x)W 2a^ ^ r^^2a^(2-e^ 

^ 2 ^ V(x) ) (e"«x - 1)2 ' ^ v{x) (e^'^x - 1)2 ' ^ 

The slow roll regime of inflation will end when e ~ 1, which correspond to the field value 

Xcnd-— ln(l + 4?) ■ (23) 



an \ V3/ 

Note that the slow roll parameter 77 is then negative, Tycnd ~ ^ ^ "73 there is a small region of negative (mass 

squared) curvature in the potential just after the end of inflation. The effective curvature of the potential will be 
negative until X* — \\i2^ which corresponds to the inflection point, given by rj^, — 0. 

During the slow-roll regime Hsr — which evaluated at = 60 e-folds is approximately given by i/go — ^■^ 

where M defines the natural inflationary energy scale of this model as well as the frequency of oscillations ([2T|) during 
reheating. At the end of inflation iJond = oi": equivalently, TJend — f-f^eo = ^■ 

The radiative corrections for a model containing a single scalar field h non-minimally coupled to gravity were 
generically calculated in Ref. The specific radiative corrections for the model under consideration, estimated in 
Ref. [ll], were recently reviewed in Ref. [2^. In what follows until the end of this section, we will summarize the 
results of this last work. For large ^ and slow varying h the main contribution comes from loops of the matter [2^[23| 
and the effective action can be calculated by a local 1 /m^-expansion in powers of the curvature and its gradients and 
the gradients of the Higgs field, obtaining [H, 

S[g^.., h]= J d^x l^~U{h) + F{h) R{g^,) - i G{h) (\/hf^ , (24) 

where the functions V{h), U{h) and G{h) are given by 



f/(/.)^>^-^^)^ + ^(Aln|+i.), (25) 
m = \m+ih^)^^^{cX.^^+D], (26) 



= 1 + 77^ (^1^752 +^) • (27) 



Here A,B_,C,D,E and F are different combinations of the Higgs, gauge and Yukawa couplings and their logarithms 
[13, HE Hal J Emd Q is the normalization scale. Following Ref. for the analysis of inflation we will just consider 
the explict form of the combination A, 



A f 



which is related with the local conformal anomaly. The factor 3 accounts for the polarizations of the gauge bosons 
and a similar factor 4 for the fermions has been taken into account. 

The new inflationary potential in the Einstein frame for a Higgs field with non-minimal coupling ^ S> 1 and mean 
value much greater than the minimum of the classical potential, is given by [23l | 



AfJ,V Uih) 



2 J F^{h) 



(29) 

h=h(x) 



The corresponding slow roll parameters are 



4M|, / h'^\'^_A(Ml A \^ AMI 
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with /ij 



6iTT^ Mi 



Using the WMAP+BAO+SN constraint ,27] at the 2a confidence level gives a value for the 

(31) 



spectral index [23| 

0.934 < ns{ko) < 0.988 
For N{ko) = 60 the value for A and for the tensor-to-scalar ratio is [1 



- 12.4 < A< 14.1 , 
0.0006 < r < 0.015, 

being the spectral scalar index running completely negligible [1^ 

- 5.6 < a X 10^ < -4.3 



(32) 
(33) 



(34) 



Let us compare this window with the one obtained from the Standard Model coupling constants at the scale M , and 
obtain the Higgs self-coupling at such scale. This differs from the analysis performed in [23| . where the coupling 
constants were evaluated at the electroweak scale instead of at the characteristic energy scale of inflation, M. At 



that scale the coupling constants for the gauge groups SU{2)l and U{\)y are roughly equal gi 



92 



0.3 and 



cos Ow = sm ow 
cos 9w, 



9\Y w l/\/2- The total anomalous scaling constant at that scale, see Eq. I|28p. in terms of 52 and 



-1 



1 



|(i + 5-C 



together with the bounds obtained from the WMAP-f-BAO-f-SN 2a c.l. constraints [23|, see Eq. 
following range for the self-coupling of the Higgs field at the scale M, 

- 0.424 < A(M) < 0.482 . 



(35) 
jive us the 

(36) 



This range can be propagated back to obtain a value at the scale through the renomalization group equations. 



A-HM) = A-i(M)-t-^log^ 



(37) 



Note that we have neglected the effect of the gauge and Yukawa couplings, since the complete solution of the renor- 
malization group equations is out of the scope of this preliminary study. If we integrate this equation and take into 
account the present observational bounds [28| , we obtain a Higgs mass in the range 



114.5 GeV < mn < 275 GeV , 



(38) 



safely within the detection range of the Large Hadron Collider at CERN. In the absence of an actual measurement of 
the Higgs self-coupling A, for the analysis of sections llVI and thereafter, we will just take different values compatible 
with the above range 



III. THE STANDARD MODEL MATTER SECTOR IN THE EINSTEIN FRAME 



The length scales are conventionally defined in such a way that elementary particle masses are the same for all 
times and in all places. This implies for instance that if under a conformal transformation the Lagrangian of a free 
particle transforms as 



mds 



C 



IP 



m ~ 



(39) 



the mass should be accordingly redefined as m = ^ to express it in the new system of units. The previous argument 
applies also for classical fields [1^. The rescaling of all fields (including the metric tensor) with an arbitrary space- 
time dependent factor fi, taken with a proper conformal weight for each field, will leave the physics unaffected. The 
physical interpretation of this symmetry is clear: It changes all dimensional quantities (lengths, masses, etc.) in every 
point of the space-time leaving their ratios unchanged. 
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In this section we will apply the previous prescription to the different sectors of the action ([3]) . Consider for instance 
the Spontaneous Symmetry Breaking sector, responsible of the masses of the intermediate gauge bosons W and Z 

SsSB D-J d-'x^^m^^W+W^- + ^mlZ^Z^^} . (40) 

In the Standard Model, the masses of the SU (2) bosons are due to a Spontaneous Symmetry Breaking mechanism, 
realized by the constant vev of the Higgs field, and therefore are constant. In our case the Higgs field evolves with 
time giving rise to variable effective masses for the gauge bosons 

mw = ^ , mz ^ — , 41) 

2 cos Ow 

with 6w the Weinberg angle defined as 9w — tan^^ (51/(72) and gi and 52 are the coupling constants corresponding 
to U{1)y and SU{2)]^ at the scale M, where the relevant physical processes during preheating will take place. As 
mentioned before, numerically this correspond to a value « f;| « 0.30, which implies sin^ 9w — cos^ 9w ~ l/\/2- 
From now on we will use these values for numerical estimations. In agreement with the above prescription for 
transforming masses and fields, the action (j40p preserves its form under the conformal transformation 



SSB 



D-J d^x^g[mlrW+W^' + im|Z^Z^} , (42) 



providing that we redefine the fields and masses with the corresponding conformal weights as 

= f7 ' ' "^^-^2- 4^ ' "'^"COS^^V^- ^^^> 

The same can be applied to the interactions between fermions and gauge bosons. Let us consider, for instance, 

Sf = Snc + Scc^ I d^x^\^W+J- + %W-J+ + -Jl-Z^JiA , (44) 



lV2 V2 ^ ^ cosOw' 

where J~ = di^j^^u^, J+ = u^j^d]^ , are the charged currents carrying the information about the couplings of the 
to the Standard Model fermions, and 



= -ULl^Ui^ - -dLl^dL TT^ULl^UL + ^-^dLT'^^L • (45) 



is the neutral current with the information of the couplings of the Z boson. In the Einstein frame the action 
preserves its form 

Sf = + SEc^j d'x^g^^W^J- + ^W-J: + ^Z,rz] , (46) 



as long as we redefine the currents as 



7^ 7± 



which is equivalent to redefine the Dirac fields 



~ d u , , 

d=7^., ^=7^- (48) 



Finally, concerning the Yukawa sector, we have, for a given family of the quark sector 



Sy D — d'^x^/ —g^ruddd + m„uu| , (49) 



where d and u stand for down- and up-type quarks respectively. The effective masses in the Jordan frame, rnf = , 
become 
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in the Einstein frame. 

On the other hand, the total decay widths, summing over all the allowed decay channels in the Standard Model 
of the and Z bosons into any pair of fermions and over all the polarizations of the gauge bosons, are given 
respectively by [13] 

where Lips denotes the Lorentz invariant phase-space factors 

Lips = ^ — — sin^ 9^ + sin"* 9^ . (52) 

The decay rates will preserve their functional form in the Einstein frame, being only changed through the conformal 
transformation of the masses, i.e. 

E _ ig^mw _ ^!^Mp ( |^|\V2_3cos30^ „ 

^^^^-^^-^l^V'^ ) - 2Lips ^^-^^ 

where 'm^r± , rhz oc (l — e""'^''*^') ^ , are the dynamical masses in the Einstein frame, see Eq. (|43p . 

IV. REHEATING IN THE STANDARD MODEL OF PARTICLE PHYSICS 

In this section we will analize the different mechanisms that could give rise to efficient reheating of the Universe, 
both perturbative and non-perturbative. The natural mechanism, given the strength of the interactions of the Higgs 
boson with the Standard Model particles, would be a perturbative reheating process right after the end of slow-roll. 
However, as we will see, perturbative reheating is not efficient enough, and non-perturbative effects must be taken 
into account. Given the shape of the potential (j20[) . very different (p)reheating mechanisms could in principle take 
place, from a tachyonic production in the region between the end of inflation and the inflection point 3l|, 3^ . [33l | , an 



instant preheating mechanism [34] and a parametric resonance effect around the minimum of the potential |35l . l36l . 131 
Therefore, it will crucial to disantangle the contribution of each mechanism and quantify their relative importance. 

In order to study the different reheating mechanisms it will be useful to have an approximate expression for the 
evolution of the inflaton. As mentioned before, we can expand potential (PO)) for the range of interest, as 

V{x) = ^M^X' + AV^(x) , (54) 

with = XMp/3£,^ the typical frequency of oscillation. The first terms of the corrections AV to the quadratic 
potential are given explicitely by 

AVix) = -^\x\' + ^x' + Oi\x\'), (55) 

with P = AMp/V6$^ and C = IXjTI^. The Klein-Gordon equation for the inflaton, 

X + 3i/x + V^'(x) = 0, (56) 

can then be written, for a power-law evolution a oc t^, as 

t^X + Sptx + t^M^[l + dM^{x)]x = , (57) 

where we have neglected the Higgs' interactions with other fields, because they are proportional to the number density 
of particles of a given species and therefore they will be negligible (see section |VT|) during the first oscillations of the 
Higgs field. The backreaction of other particles into the dynamics of Higgs will only be relevant once their occupations 
numbers have grown sufficiently. The non-linear terms of the Higgs' self-interaction, described by 

6M^^--(3\x\+Cx'+0{x'), (58) 
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will be also negligible from the very beginning of reheating, |(5M^(x)| ^ 1, as we will justify a posteriori. Thus, 
neglecting such a term in the effective equation of Xi the general solution can be expressed as 

Xit) = [AJ+,{Mt) + BJ^,{Mt)] , (59) 

with A and B constants depending on the initial conditions (end of inflation), and J±^{x) Bessel functions of order 
±i/, with V = (3p — l)/2. For a reasonable power index, p > 1/3 - for matter p = 2/3, while for radiation p = 1/2 - 
the second term in the right-hand side of Eq. ([55]) diverges in the limit Mt and therefore should be discarded on 
physical grounds. The physical solution is then simply given by 

(3p-l) 

X{t) = A (Mt) — J(3p-i) {Mt) , (60) 

2 

which making use of the large argument expansion {Mt 3> 1) of fractional Bessel functions (38j . can be approximated 
by a cosinusoidal function 

X{t) - A^{Mtr"-^ cos {Mt - (3p/2)(7r/2)) . (61) 

The normalization constant A can be fixed if we consider that the oscillatory behaviour starts just at the end of 
inflation, i.e. x{t = 0) = Xond = Q;~^Mp log(l + 2/\/3) ((23| . In this case, 

A = Xe„d2^(3p-i) Q(3p_l)^,^ (62) 
where we have made use of the limit of the Bessel functions when Mt <C 1 

Jin„ u{Mt) w . (63) 

2(3p-i)V J 25(3p-i) (i(3p-l))! ^ ^ 

The energy and pressure densities associated to the general solution (|60|) are given, after averaging over several 
oscillations, by 

Px + ^^^'^') ~ lM^X^[{cos\Mt - 3^p/4)) + {sm^{Mt - Sirp/i))] = ^M^X^ , (64) 

Px ^ (^X' - « hl^X^[{cos\Mt - 3^p/4)) - (sin2(Aft - 3^p/4))] = , (65) 

with X{t) (X {Mt)~^ . Since the averaged pressure is negligible p^ ~ 0, then a{t) cx t^ with p w 2/3. Using this fact, 
the physical solution is finally expressed as 

x(t) = ^sin(Aft). (66) 

Rewriting the previous equation in terms of the number of times the inflaton crosses zero, j = (Mt)/7r, or equivalently 
in terms of the number of oscillations N = j/2, then 

X{t) « 1^ sin(27r7V) = ^ sin(7rj) ^ X{j) sin(^j) , (67) 

Therefore, the Higgs condensate oscillates with a decreasing amplitude X{j) cx l/j. We can obtain an upper bound, 
aK\x\ < 0.122/A^, on the amplitude of the Higgs field after N oscillations, which in terms of the correction of SM"^ to 
1, see Eq. implies ISM"^] < 0.122,0.0615 or 0.0244, after the first iV = 1,2 and 5 oscillations, respectively. Thus, 
from the very beginning, the effective potential of the Higgs field tends very rapidly to that of a harmonic oscillator, 
which justifies a posteriori the approximation |(5M^| <C 1 used in the derivation of Eqs. ([59|l and ((67|l . 

Note that if we neglect the presence of other fields and consider the Higgs-condensate as a free field only damped 
by the expansion rate, then we can easily estimate the number of semi-oscillations before the amplitude of the field 
becomes smaller than the transition value xt ~ -^p/Cj defined in section |TT1 For \x\ < Xt, the Higgs potential 
will not be anymore approximated by Eq. Ij54p . but rather by a quartic form (A/4)x^. This will happen when 
X{jt)/Xt ~ ^i^jf^ < 1, which impHes 

.>,..iMi±2M.<,ao'), (6S) 

an 

Therefore, if before such a moment, the Higgs field has not yet transferred efficiently its energy into other fields, then 
after jt semi-oscillations, the transition in the behaviour of the potential will also imply a change in the expansion 
rate, from a matter-like, characteristic of quadratic potentials, to radiation-like, characteristic of quartic potentials. 
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A. Perturbative Decay of the Higgs Field 



A natural reheating mechanism, given the strength of the interactions of the Higgs boson with the Standard Model 
particles, would be a perturbative decay process of the inflaton quanta into the Standard Model particles right after 
the end of slow-roll. As we saw, soon after the end of inflation the effective potential for the Higgs field can be 
approximated by a simple quadratic potential (|2ip . In this approximation, the masses of fermions and gauge bosons, 
see Eqs. (|43)) and (|50|) . are simply given by 

/^_a|xl V^\f ( "1^1 \^''^ nr 92 ( a\x\ \^''^ ,«n^ 

In order to have a perturbative decay two conditions must be fulfilled : 

1) There should be enough phase-space in the final states for the Higgs field to decay, i.e. M > 2mf,mA, which 
will only happen when the amplitude of the Higgs field becomes smaller than a certain critial value Xc- In particular, 
for a decay into gauge bosons and/or fermions, in the light of Eq. (|69p . one needs 

1 

where g = 32,52/ cosdw for the W and Z bosons, respectively, and g = yj for the fermions. When compared to the 
initial amplitude (l23l) of the Higgs field at the end of inflation, Xc/Xend — l/(3g^0, see that this critical value is 
much smaller than Xond for the gauge bosons and the top quark; of the same order for the bottom and charm quarks, 
and even greater for the rest of the quarks and the SM leptons. 

2) The Higgs decay rate F ^ ^^'^ to be greater than the rate of expansion = w ^(bj)^ (^^5^) ' 
where we have used Eq. (|64|) . Such a condition, T > H, can be translated into the following inequality 

J>J.^^^^^, (71) 

which defines the critical number of semi-oscillations required for this second condition to be true (again, 
9 — 92, 92/ cos^u) and y/, for W, Z bosons and fermions, respectively). 



X>Xc^-\hM, (70) 



The critical amplitude ((70|) below which the Higgs is allowed to decay into gauge bosons is of order Xc ^ O.IM. 
As mentioned before, this amplitude is much smaller than that of the Higgs at the end of inflation, Xond ~ 0.94Mj, w 
8 X lO^M. Therefore, the Higgs condensate would need to oscillate ^ 10® times before being able to decay through 
this channel. The same applies to the top quark. In the case of other fermions, due to the wide range of the Yukawa 
couplings, several situations can take place. For instance, the decay channel into bottom and charm quarks is opened 
only after a few oscillations of the Higgs, while for the rest of quarks and leptons, the decay-channel has sufficient 
phase space from the very end of inflation. In general, the smaller the Yukawa coupling of a given fermion species 
to the Higgs, the less oscillations the Higgs will go through before there is enough phase-space for it to decay into 
such fermion species. Notice however that the smallness of the Yukawa coupling implies also a smaller decay rate. 
Consider for instance the decay of the Higgs into electrons, whose Yukawa coupling is of order ye « 10~®. From the 
very end of inflation, see Eq. (|70p . there is phase-space in this channel for the Higgs to decay into. However, it is 
precisely the smallness of the electron's Yukawa coupling that allows the decay to be possible, which prevents the 
condition 2) to be fulfilled. The decay width is much smaller than the Hubble rate for a huge number of oscillations. 
Looking at Eq. ([7T|) . we realize that the Higgs condensate should oscillate j ~ IQ-'^^ times before the decay rate into 
electrons overtakes the Hubble rate. 

One can check that the previous conclusions also hold for the rest of fermions of the Standard Model. When there 
is phase-space for the Higgs to decay into a given species, the decay rate does not catch up with the expansion rate 
and, viceversa, if the decay rate of a given species overtakes the expansion rate, there is no phase-space for the decay 
to happen^. Therefore, during a large number of oscillations, the Higgs field is not allowed to decay perturbatively 



^ Note that the condition I I70I I (which prevents Higgs decay into gauge bosons and top quarks) assumes an average amplitude over a single 
Higgs oscillation, while smaller values are attained around the minimum of the potential when X{t) < Xc- However, when this happens 
the Higgs field is well inside the non-adiabatic range |x| < Xa Il80|l . in which the very concept of particle (gauge bosons and top quark) 
is not properly defined, see section [IV CI 
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in any of the Standard Model fields. Moreover, before any of those decay channels is opened, many other interesting 
(non-perturbative effects) will take place, as we will describe in detail in the next sections. 



B. Tachyonic preheating and Non-adiabatic particle production at the inflection point 

As we pointed out in section |TT1 the effective square mass of the Higgs field x is negative just after the end of 
inflation and will be so till the inflection point. When this happens spinoidal instability takes place [3l|, HI, HI] and 
long wavelengths quantum fluctuations Xk, with momenta k < m-^, grow exponentially. The width of the tachyonic 
band will be limited in our case by the point of maximum particle production, the end of inflation. At this point the 
effective mass TOv takes a value 



"^^(Xcnd) = 



30 



(72) 



which corresponds to a maximum momentum for the tachyonic band kmax — 0.2M . This comes from vacuum quantum 

exp(Mi^l/30- {k/Mf), which grow exponentially. 



fluctuations, Xk{t) cc exp^it^k^ + M^) 

However, since the inflaton is fast rolling down the potential towards the positive curvature region, the duration 
of the tachyonic preheating stage is so short that the occupation numbers of those modes in the band do not 
grow significantly and the effect can be neglected. In particular, the time interval from the end of inflation till 
the inflection point is just MAt « 0.5 and therefore, even for the fastest growing mode. A; = 0, its growth is only 
~ g0-5/\/30 ^ 1.09. This is a negligible effect and thus, one can still consider an initial spectrum of quantum vacuum 
fluctuations even at the inflection point. For simplicity, all our analitycal estimations have been done ignoring 
this period of tachyonic instability, taking as initial conditions at the end of inflation, the amplitude of the Higgs 
condensate Xend ~ log(l + 2/V3) (|23p and quantum vacuum fluctuations. 



Another physical effect before the Higgs condensate reaches the bottom of the potential for the flrst time, will be 
the particle production in the inflection point, due to the violation of the adiabaticity condition. 



\^k\ » K\ 

where the frequency of oscillation of the fluctuations is uj^ it) = 

adiabaticity condition as uJk^k ^ uj"^, we find that only those modes within the band 



(73) 

^2 ^ v"{x)- Differenciating this and rewriting the 



k' 



V {x)x 



(74) 



are amplified. At the end of infiation H = —H^ which implies x ~ ~^^^^(x)- Extrapolating the previous formula to 
the inflection point (ip) we get 



4k 



(75) 



which indeed seems to be a very good approximation if we compare it with the result of a numerical solution beyond 
slow-roll. Inserting ([75| into Eq. ([74|) we get 



K « 



V (x)Xi, 



^60 



(76) 



which corresponds to a maximum excited wave number at the inflection point given by kip < 0.4Af . Again, here the 
time of production is so brief that the occupation number of modes within the band is not signiflcantly enhanced. 
We will have to wait until the next stage of consecutive oscillations of the Higgs field around zero for a significant 
production of particles. 



C. Instant Preheating 



During each oscillation of the Higgs field x, the rest of the quantum fields that couple to it will oscillate many times. 
Consider for instance the interaction of the Higgs field with the Z bosons around the minimum of the potential. In 
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this region the associated action (I42p can be approximated by a trihnear interaction where the masses of the W and 
Z bosons (|43|) are given by 



9w 



^W-^^\X\, ^^" 4ecos^g, (^^^ 



which are much greater than the inflaton mass AI for the main part of the oscillation of x. As a result, the typical 
frequency of oscillation of the gauge boson is much higher than the one of the Higgs field x- This implies that during 
most of the time the effective masses of the intermediate boson are changing adiabatically and an adiabatic invariant 
can be defined: the number of particles. However, for values of x very close to zero, the adiabaticity conditions 



mw 



''W 



< 1 



mz 



ml 



« 1 , (78) 



are violated. In such a case, there will be an inequivalence between the vacua before and after the passage of x 
through the minimum of the potential, which can be interpreted as particle production [35l . [36j . In terms of the field 
X, the violation of ([75)1 corresponds to the region —Xa ^ X ^ Xa , 



1/3 



\agmpj 

where, from now on, g = g2, 52/ cosOw for the W ot Z bosons. Only outside this region, the notion of particle makes 
sense and an adiabatic invariants can be defined. Taking into account that and approximating the velocity of the field 
around zero as xU) ~ M — MX{j), see Eq. (|67| . the general expressions ([79|) can be approximated as 



Xaj^g-^TT-^Mp J V3Clog(l + 2/V3)/ gs 

Note that the previous regions are indeed very narrow compared to the amplitude of the oscillating Higgs, Xa ~ 
\Q~'^i^/^X{j). Therefore, the particle production that takes place in that region happens within a very short period 
of time as compared to the inflatons' oscillation period T — 27r/M, 

Ata(j) ~ % ^ 10-2 //3 « T , (81) 

Ixl 

Notice indeed that different values of A do not change appreciably the above conclusions about the smallness of the 
non-adiabatic regions. Given the weak dependence of Ai cx j^^'^, many semi-oscillations {'-^ 10'^) will pass before the 
fraction of time spent in the non-adiabatic zone will increase from a 1% to a 10%, as compared with the period of 
oscillations. This holds also independently of the species, W ot Z bosons. 

Moreover, despite the smallness of Xo as compared to the amplitude X{j), it is important to note that the field range 
corresponding to the region of non-adiabaticity still is much greater than those critical regions defined in section [Hi 
In particular, let us recall that there is a field value, Xt ~ ^p/£., below which there is a transition of the effective 
potential from a quadratic to quartic behaviour. However, this is well inside the region of non-adiabaticity, Xt ^ Xa, 
as we emphasized before. Moreover, there is also an interval of Higgs field values, |x| < Xc for which the Higgs 
perturbative decay into W, Z and top quarks can occur [see Eq. (|70|) in section flVAj . which nevertheless is also much 
smaller than the non-adiabaticity interval, Xc ^ Xa- 

We will now discuss the non-perturbative creation of particles in the non-adiabatic region. This production is 
formally equivalent to the quantum mechanical problem of a particle scattering in a periodic potential. In the case 
under consideration the equations of motion for the fluctuations of each gauge field with a given polarization will be 
given by VF^' -|- (fc^/a^ + rn^)Wk = 0, and the corresponding one for the Z-fluctuations. Expanding Eq. (|57| around 
the j-th zero at time tj — nj, the evolution equation of the fluctuations can be approximated as 



/A^ .g|M,Xcnd|sin(M(^-t,0)h 
W 47rj^cos2^ ) 

Notice that around the zeros of the inflaton the sinusoidal behaviour | sin(A/(t — tj))\ can be very well approximated 
by I sin(M(t — tj))\ « M\t ~tj\ = r, which allows us to rewrite Eqs. ([5^ and ([55)) as a Schrodinger-like equation like 

MM/fe^i^'VFfe, -4'--M^fc-^'^fc, (84) 
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FIG. 2: Spectral distributions (|88p for the gauge bosons created in a single zero crossing through the first term of Eq. (|86l) . 
calculated after j = 1, 2, 5 and 10 oscilations (from left to right). The horizontal axis represents x = k/Mq^''^ , so x = 1 is the 
typical width of the band of momenta of particles created at the first scattering. For later times, the distributions broaden out 
to greater momenta, since the argument of Eq. (|88p . Xj behaves as oc The typical momenta of the distribution agree 

with the one calculated in section HV CI 



where primes denote derivatives with respect to the rescaled time t = Mt, K is the rescaled momentum K = and 

q^^COse^q, = ^-^^ — j = , (85) 

are the usual resonance parameters fSG^. Each time the inflaton crosses zero can be interpreted therefore as the 
quantum mechanical scattering problem of a particle crossing an inverted triangular potential. Let T and R = \ — T , 
for either W or be the transmission and reflection probabilities for a single scattering in this periodic triangular 
barrier. The number of particles just after the j-th scattering, n}.{j~^), in terms of the previous number of particles 
nk{j~) just before that scattering, can be written as [ssl. [3^ 



in = iT-\j) - 1) + {2T-\j) - l)n,(r) + 2cose,^T^\j) {T^\j) - l)^nk{r) K(r) + 1) , (86) 

where 6j are some accummulated phases at each scattering, that we will discuss later on in section |Vl since they will 
not play any role in the following discussion of this section. The inverse of the transmission probability for the j-th 
scattering can be expressed as [30| 

Ti:\j) = l + n' [Ai Ai' (-x,2) + Bi Bi' i-.'^)]' , (87) 

with Xj = K/{q/ jY^^ and Ai(z),Bi(z) the Airy functions. Note that we have used the Wronskian normalization, 
Ai(z)Bi'(z) - Bi(z)Ai'(z) = n-\ 

Consider the situation, nk{j~) <C 1, which is certainly true in the first scattering j = 1, or can happen for > 1 if 
the previously produced gauge bosons have fully decayed into fermions. In such a case, 

Anfe(j + ) « T-\j) 1 , (88) 

where we have retained only the first term of Eq. ([55]) . This corresponds to the spontaneous particle creation of W 
and Z bosons each time the Higgs crosses zero and, therefore, tells us about the number of particles of these species 
that are created in each zero-crossing. The momenta distribution is shown in Fig. O In particular, the total number 
of produced particles of a given species with a given polarization, just after exiting the non-adiabatic region around 
the j-th zero-crossing, can be obtained as 

1 

= / dk e T^' U) - 1 =^IM\ (89) 







2j 



with 1 — [Ai(— a;^)Ai'(— a;^) + Bi(— a:^)Bi'(— x^)] ^ x^rfx ~ 0.0046 and q the resonant parameters given by Eq. ((85|) . 
Thus, the only difference between the number of W and Z bosons produced is simply encoded in the different resonance 
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parameter, qw cx: and qz cx; g2 / cos^ 9w , respectively. Notice that the effect of the non-perturbative production of 
these particles is proportional to the coupling square. If the couplings of the Higgs field to the gauge bosons were not 
so large (52 0.5, cos~^ Ow ~ 1-4), then their produccion would be very supressed. Strictly speaking the previous 
analysis is just valid for gauge bosons. The production of fermions through this mechanism is different and more 
involved than for bosons. Nevertheless, if the effect is, as expected, proportional to the Yukawa coupling squared [i^ 
then only the top quark production would be non-negligible. 

After the passage through the minimum of the effective potential the number of W and Z particles remains (almost) 
constant while their masses grow when the field x increases. The W and Z bosons tend to decay into fermions in 
a time At ^ {^wz)J^^ where z given by Eq. ([53]), while represents a time average between the j- and 
the (j + l)-th scatterings. Given the time-dependence of the x field ((57)) . the typical time of decay turns out to be 
At ~ 0.64 j^/^ for the Z bosons and a bit bigger, At ~ 1.55 j^^^ M~^, for the W bosons, as was expected. This 
implies that in a semiperiod, T/2 = ttM~^, the non-perturbatively produced gauge bosons at the j-th scattering, 
decay significantly before the next scattering takes place, at least for the first scatterings. However, as the amplitude 
of the Higgs field decreases with time due to the expansion of the Universe, the probability of decay of the gauge 
bosons, sec Eq. ([55)) . becomes smaller and smaller as time goes by. This explains the j^/^ behaviour, which essentially 
means that after a certain number of oscillations, the number of produced fermions through the perturbative decays 
per semi-oscillation, will become eventually negligible. 

In the Jordan frame the Standard Model presents its usual form and the fermions produced in the decay of the Z 
and W bosons are mainly relativistic. Since both momenta and masses transform in the same way under a change of 
conformal frame, if a gauge boson is allowed to decay into a pair of fermions in the Jordan frame, it will also be able 
to decay in the Einstein frame. Therefore, the relation between the typical momenta and masses of those fermions 
(F) and gauge bosons (W,Z) in the conformally transformed frame, is simply given by 

2Ckp + ml) ^ klr + ml, , 2(fc| -t- m^) = fc| + m| . (90) 

In terms of the field Xi the previous equations can be rewritten as 

(91) 

- 1 , (92) 

for the W and Z fields respectively. Note that the relativistic or non-relativistic nature of a given particle is something 
intrinsic to the particle and should not depend on the conformal frame. As expected the transitions Z tt, W tb 
are not allowed in the Einstein frame. For the rest of quarks kp 3> mjp, which implies that all the fermions produced 
in the decay of the W and Z bosons are clearly relativistic, as happened in the Jordan frame. The total number 
density of gauge bosons n{j~^) present just after the j-th crossing will decay exponentially fast until the next crossing, 
due to the perturbative decay into fermions. Therefore the total number density just previous to the (j + l)-th zero 
crossing, n((j -|- 1)~), is given by 

niii + 1)") = n{j+)e~^'^*' = 7i(j + )e-<r>^ i . (93) 
The number of fermions produced between those two scatterings, Anpij), is simply given by 

AnF(j) = 2x3x L2(j+)(l-e-<r^>^*) + 2nH^(j+)(l-e-<^"->^*)j , (94) 

where the factor 2x3 takes into account that each gauge boson can have one out of three polarizations and decay 
into two fermions, while the extra factor 2 in front nw, accounts both for the and W~ decays. The averaged 
value of the decay widths in the previous expressions can be estimated, see Eqs. ()53p . as 





1 






k'p 


1 


mj. 






/r \ 3 cos'^ Ow ,^ , _ p. .„„^ 
{Tw^aiOj = ^j^.pg {Tz-^aii)j = -jrF{j) , (96) 



where T — 2Tr/M is the typical oscillation period and we have defined 




0.3423^. (97) 
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Note that the last approximated equahty is simply a (good) fit to F(j) for all j. The constants JZtJw are just 
numerical factors depending of the parameters of the model and the decaying species, 

iz = ( ^ Lips - 14.23A-^ , . ^|^7z - 5.91A-^ (98) 

\cos Ow J 16A^/-^ zLips 

Using the notation, Epzij) and Ep^{j) for the mean energy of the fcrmions produced between tj and tj+i, from 
the decay oi Z or W bosons, respectively, then we find 

Ep, (j) ^ ^ ^kl + ml'j « (kp), « 1 (mz), « ^^^/^ZsBw ^^^^ ' ^^^^ 

Ef^U) = l^^kl + ml^ « {kF)j « \{mw), ~ ^ F{j) Mp , (100) 

where we have used the fact that the produced fermions are relativistic, see Eqs (PT|) and ([5^ . while the gauge 
bosons are non-relativistic, see Eq. p30p in section IVl 



Let us work now under the following hypothesis: we will consider that the perturbative decay of the gauge bosons 
into fermions is sufficiently effective, such that the gauge bosons do not accumulate significantly. This amounts to 
neglect initially a potential effect of parametric resonance. This is of course a rough approximation, which is only 
valid for the first oscillations, where e"'*'^*^^^ <C 1. Numerically, after j = 1,2,10,15 and 20 zero-crossings, the 
99.5%, 98.5%, 94.2%, 87.4%, 81.9%, 77.4% respectively, of the produced Z particles have decayed into fermions (and 
a similar though smaller fraction of the W bosons) . This implies that there will be always a remnant of the gauge 
bosons produced at each scattering, that will not decay in one semi-period of the inflaton's oscillation. Let us neglect 
this for the time being, therefore ignoring the possibility of having parametric resonance, and estimate the energy 
transferred simply through the perturbative decay into fcrmions, during the first oscillations. 

In particular, the energy density of those fermions produced after the first scattering, averaged over the first semi- 
oscilation between Mt = it and Mt — 2tt, will be 

Ap;^(l) ^ 6 [Anz{l)Ep,{l) + 2Anw{l)EF„{l)] = e (Ul^l^^^) F{1) , (101) 



2 "-^^TT^ 



with 



£ = 



_ {2 + cos-^ 0w)lg; 



3 X IQ-^X-^/^ . (102) 



The energy density of the inflaton, evaluated at the maximum amplitude of the first semi-oscillation, is given (j64p by 

Px(l) « ^M^xLd (1^) • (103) 
Therefore, the ratio between the energy density of the fermions and of the inflaton, at that moment, Aft sa 1.57r, is 



e{l)^^^^^2.10-^X^y\ (104) 
Px(l) (2/3)^ 

which means that initially, for e.g. for A = 0.4, only a ~ 0.004% of the inflaton's energy has been transferred to the 
fermions. Thus, the so called Instant Preheating mechanism ([sj]) results frustrated here, because in order to make 
it work efficiently, the couplings of the theory must be really fine-tuned, in such a way that a significant fraction of 
the energy of the inflaton was transferred (in the first semi-oscillation) to the decay products of the bosons to which 
the inflaton is coupled. Moreover, in the instant preheating scenario, the produced fermions must be non-relativistic 
while the effective behaviour of the background inflaton should be effectively mimicing that of relativistic matter (like 
e.g. in Ax^ models). Only in this case it would be guaranteed that the remnant energy of the inflaton would decay 
faster that that of the fermions, thanks to the extra supression factor 1/a due to the expansion of the Universe. If the 
inflaton would effectively behave as non-relativistic matter and the produced fermions were relativistic, the energy of 
the inflaton could again overtake very soon that of the fermions, because now the fermion's energy would decrease 
faster than that of the background. That is, precisely, the situation we have in the scenario under discussion. Even 
if we had found that e(l) ~ 0{1), the relativistic nature of the fermions and the non-relativistic effective behaviour 
of the Higgs oscillations, would have prevented the Universe to instantaneously reheat at that point. 
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One could hope that after a certain number of oscillations, let us say jp, that ratio would grow up to a value 
e{jp) ^ 0{1). The succesively produced fermions, generated each semi-oscillation through the perturbative decay of 
the (non-perturbatively produced) W and Z bosons, could perhaps accumulate a sufficiently amniount of energy that 
could finally equal that of the Higgs condensate. This seems not totally unresonable because the total energy stored 
in the Higgs decreases with the expansion of the Universe as oc 1/j^, see Eq. ([M)) . such that the total amount of 
energy that we would require to transfer to the fermions would be less and less. Moreover, the number of fermions 
would only increase as time goes on, so one keeps adding energy each semi-oscillation in the form of new produced 
fermions. Thus, these two effects would contribute to the increment of the ratio of the energy between the fermions 
and the Higgs. On the other hand, the relativistic nature of the fermions and the decrease of their production rate 
with the expansion, as F cx 1/Vj, would tend to decrease such a ratio. Therefore, one must put together all these 
competing effects in order to obtain the evolution in time of the energy transferred from the Higgs to the fermions. 
To do this, we will assume, both for simplicity as well as for trying to make this mechanism more efficient, that since 
the gauge bosons do not accumulate significantly, only the first term of Eq. (jSG]) should be considered, such that W 
and Z bosons are only produced through spontaneous creation at each zero-crossing. 

In this case the averaged energy density of the fermions produced between tj and tj+i will be given by 



where e is given by Eq. (|102p and we have defined 



(105) 



T(j) ^ ( 1 - U + ^co^ >^we 1 \ 

^' \ (l-|-2cOs3 6'vy) J ^ ' 

Then the ratio between the energy of the fermions to the Higgs condensate at the j-th zero crossing, finally reads 

e(.) ^ ^ - t {^-e GU)il±^ , (107) 

PxU ) J \J / J 



where 



G(J)^^^^(^)T(^)(^)'^^ (108) 

Note that the strength of the effect, i.e. the amplitude of e{j), is modulated by the gauge couplings through e oc 
so, even in this case in which the SM gauge couplings of the Higgs to the vector bosons are quite big (5I ~ 0.3), 
that does not help to transfer sufficient energy initially. As mentioned before, if we could apply this formalism to 
the production of fermions at each Higgs' zero crossing, by substituying the gauge couplings with the Yukawa ones, 
we would obtain even a more ridiculous production of particles (except perhaps for the top quarks) . Of course, the 
question of fermionic preheating at each zero crossing deserves more investigation and we will address it in a future 
publication. 

The numerical values of the ratio e{j) after e.g. j = 1,2,5,10,15 and 20 semi-oscillations, for e.g. for A — 0.4, 
are respectively e(j)[xl0^] 3.90,5.97,11.82,37.26,65.04 and 97.59. For different values of A, these numbers do 
not change significantly. Thus, we see that the transferred energy from the Higgs field to the fermions through the 
gauge bosons is generically a very slowly growing function. After 20 crossings the transferred energy is still only 
~ 0.03% of the Higgs energy at that time. Therefore, we clearly see that this succesive Instant Preheating mechanism 
is not efficient enough as to rapidly reheat the Universe. If we consider that the former formalism (|107p is valid up 
to an arbitrary number of oscillations, then we can estimate the number jp of semi-oscllations required to achieve 
e(j) ~ 0{1). Equating e{j) to 1 in Eq. ()107p . we obtain jp ~ 0(10**). However, much earlier than that, parametric 
resonance effects should be considered, see section IVl 

In other words, notice that we have neglected the presence of those Z and W bosons that did not decay into 
fermions in each semi-oscillation. The occupation number of the bosons produced at the bottom of the potential is 
not simply generated by the first term of Eq. but rather by the rest of the terms in Eq. ([5^ . which indeed 

give rise to the phenomena of resonant production of bosons. Taking this into account will have very interesting 
consequences. The number density of the bosonic species will grow exponentially fast and thus will also transfer 
energy into the fermions exponentially rapidly. We must therefore develop a mixed formalism that takes into account 
the two competing effects: that of parametric resonant production of bosons versus the effect of their perturbative 
decay into fermions. It is crucial to note that while the perturbative decay does not transfer enough energy (as we 
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have just seen), the fact that those bosons disappear wiU have very important consequences for the development of 
the resonant effect. In particular, the resonance will not become effective from the beginning of the oscillations of 
the inflaton right after inflation, as usually assumed, but only after the inflaton has already performed a significant 
number of oscillations. 



COMBINED PREHEATING: MIXED PARAMETRIC RESONANCE AND PERTURBATIVE 

DECAYS 



Let us now analyze how the occupation number of the Z and W bosons grow if we consider the effect of all the 
terms in Eq. (|86p . The production of gauge bosons will also occur, as before, in a very short interval of time (|8ip 
when the Higgs condensate crosses around zero, violating then the adiabaticity conditions |x| < Xa- In the 

large occupation limit rt^: ^ 1, the first term in Eq. ((86)) can be neglected and therefore the spectra number density 
of the produced gauge bosons just after the j-th scattering is given by 



(109) 



which indicates the spectral number density nk{j^) just after the j-th scattering, in terms of the spectral number 
density nk{j~) just before such scattering. Since the interval between successive scatterings is AfAi = tt, we can 
define naturally a growth (Floquet) index /ife(j) as [ssl [36t 



nfc(j + ) « nfc(j-)e2'^'=(^')^^^* = nfe(r )e'"^^^^'^ • 



Comparing formulas, we obtain 



/^fcW«^log((2T,7\j)-l) 



2 cos 9 



(110) 



(111) 



The Oj are some accummulated phases at the j-th scattering, which can indeed play a very important role, since they 
can enhance (cos 9j < 0) or decrease (cos 9j > 0) the effect of production of particles at each scattering. 

Depending on the phases, we can consider the following cases: The typical behaviour of the Floquet index, for 
cos 9 = 0, 



„(typ) 



^log(2T--l), 



(112) 



the maximum index, aehieved for eos^ = — 1, given by 



(max) 



1 



■log 



T7 



and the average index over an oscillation, obtained as 



av) 



. i- £%.(.) ...^ log (2T-) 



(113) 



(114) 



All these possibilities are shown in Fig. [31 as a function of x = K / {q/ j)'^/'^ , where q are the resonant parameters ([55)1 
for the Z and W bosons, while Xj is the natural argument of the transmission probability scattering functions ([87]). 



As explained in Ref. [36|, when A9j = 9j+i — 9j 3> tt, the effect of resonance will be chaotic, being then the phases 
essentially random at each scattering. For instance, using the effective frequencies of the fluctuations ([82]) of the W 
field, these phases can be estimated, for the relevant range of momenta, as follows 



2VA 



F(j) ~ O(10-2)j 



-1/2 



(115) 



where -F'(j) was defined in Eq. (j97p and we have neglected K versus in the second equality, since, as will be 
justified later p30p . the produced bosons are non-relativistic. Comparing the above formula with tt, we see that the 
end of the stochastic behaviour will occur after ~ (2 — 5) x 10"^ zero crossings, depending on A. For the case of the 
Z boson the previous estimation of the end of the stochastic resonance is modified by a factor (cos6'w)^^ ~ 
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being thus the result essentiaUy unaffected. Therefore, since for the first thousand of oscillations of the Higgs, the 
accummulated phases of the fluctuations of the gauge bosons will be chaotic, we will average out the phases and work 
with ^i'^^^ 

On the other hand, the pcrturbative decay of the produced vector bosons occurs precisely just between two successive 
Higgs zero-crossings, n{{j + 1)^) = exp(— 7 where F{j) is given by Eq. ([97]) and 7 — 'yzilw, see Eq. 

([98|) . Taking into account Eq. ()110|1 and Eq. ()111|) we can express the number of gauge bosons just after the (j + l)-th 
scattering in terms of the number just after the previous one 

+ 1)+) = M{3 + l)-)e'"^''(^+i) = nfc(j+)e-^^(^')e2"^'=(^'+i) , (116) 

Applied recursively, this formula allows us to obtain the occupation number for each species and polarization, just 
after the [j + l)-th scattering in terms of the initial abundances nfc(l"'"), 

j 3 

+ ='^fe(l+)cxp [-7^^(0] exp [2^^Aife(* + l)] ■ (117) 

i=l 1=1 

The initial abundances are, of course, only generated through Eq. (|88p . and are given by 

nfe(l+) = T^\l) - 1 = n^Ci'i-xi) , (118) 

where we have defined the function 

Ciixj) = Ai(-a;2)Ai'(-x2) + Bi(-a;2)Bi'(-a;2) , (119) 



Mq^/^aj 



(120) 



Again, we used Xj in light of Eq. (j84p . as the natural argument of the expression of the transmission probability (j87p . 
Note that here, the species are only distinguished through the resonance parameters in Eq. (j85p . Normalizing the 
scale factor at the first zero crossing as ai = 1, then we can simply write the evolution of the scale factor as aj = . 
Thus, the behaviour of Xj with the number of zero crossings goes as oc Then, we can define a typical momentum 

of the problem, A:*(j), related in a very simple way to the resonance parameters qz^ qw (|85p . as 

^ = (121) 



' jl/3fc,(l) 

where 

fc.(l) ^ q'^'M ^ (^^^^f^) ' M , (122) 

with g = §2 and 52/ cos 9w for W and Z bosons, respectively. Since fc*(j) is the natural scale for the momenta of 
the problem, its order of magnitude should coincide simply with the one obtained via the Heisenberg uncertainty 
principle, see Eq. ([8T|l . as it is indeed the case since 

k{j) ^ a.iAta)-' ^ « K{3) . (123) 

Notice that the typical momenta range will be red-shifted because of the expansion of the Universe and, even the 
comoving typical moment fc*, is not a static quantity but rather depends on j. 

Let us obtain now the total number density of created particles. Just after the j-th scattering, this will be given by 

n{j+>2) - -L^ei-''^^^^ /dfcfcV(l+)e{'"^-^'"=^*)} 
27r^a| J 

= ^3^g{-7ECl n^} f duu^Gi\~u') f\{27r^Ci'i~u'/f/'') + 1) , (124) 

where fc*(l)'^ cx g"^ ()122|) should be evaluated with 5 = 52 or (72/ cosOw, and 7 = jw or jz, respectively for W or Z 
bosons. This formula encodes the usual resonant behaviour discovered in the 90's, see Refs. [H,!!^, in which it was 
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FIG. 3: Left: The Floquet index for a given polarization of the W and Z bosons as a function of the variable Xj — k/kt{j). 
Here we show the maximum (continuous red), the average (short dashed green) and the typical (long dashed blue) indices. 
Right: The initial spectral distribution nfe(l^) (lower blue curve) and the Gaussian approximation n(j^ > 2) (|127p for different 
j's greater than 2 (rest of the curves), describing the resonant behaviour. The approximation is so good that it is hard to 
distinguish it from the real curve, presenting small deviations just on the tail. The horizontal axis is x = fc/fc«(l) and the 
curves correspond to different j's. It is clearly distinguishable the fact that only the range a; < 1 (fc < fc*(l)) is filtered and 
therefore excited through parametric resonance, no matter if j ^ 2. 



implicitely assumed that the produced bosons didn't decay between succesive inflaton zero-crossings. However, as we 
saw in section llV C| the bosons produced each time the Higgs condensate crosses zero, significantly decay before the 
next scattering. Therefore, we had to correct our formulas for this effect. Fortunately, this was easily done, since 
the resonant growth occurs in a step-like form, instantaneously (within a time A^q, see Eq. (|8T|)) when the Higgs 
condensate crosses around zero, while the perturbative decay of the produced vector bosons occurs during the time 
just between two successive Higgs zero-crossings. Thus, the occupation number just before the (j -I- l)-th scattering, 
in terms of the occupation number just after the j-th scattering, has been corrected by the factor exp{— 7 X]i=i 
which accounts for the accummulated effect of the perturbative decays up to the j-th scattering. 

The combined effect of the non-perturbative parametric resonant at the non-adiabatic regions at the bottom of the 
potential, together with the perturbative decay along the adiabatic zone during the rest of the semi-oscillation, give rise 
to a new phenemenology, as we will inmediately see. Therefore, to emphasize the difference from the usual parametric 
resonance or instant preheating- like mechanisms, we will call these effect Combined Preheating. Expanding the 
combination of Airy functions (|119p . it is possible to write 

Ci2(x,)«Ce-^"'/^"'', and i27:'Ci\-uyf /') + 1) ^ Ae-""^" (125) 
where u = j^^^xj and 

r (2/3)^ n-JL&l ._o,o^V R 1 (16/3)2^^ 1 

r2(i/3)r2(2/3)' 32/3r(i/3) ^-^ + ^^^' 3y^r^{i/3)r{2/3) A' ^ ^' 

Substituying Eqs. (fT25)) and (fT26l) in Eq. p24)) . then we obtain 

n(j+>2) « M'''"^''l~'^^''*^^h^-'C J duu^e-^"^e-^g^^^'-^^^)"^ 

= Af3f::!l|^^.-C^ [D + B±ry^Y^' . (127) 
where we have used aj — j"^^^, performed the resulting gaussian integral and defined 

(128) 

1=1 

for simplicity. Notice that the resonant behaviour is now encoded in the factor A^^^, which, for sufficiently great 
j, will finally overtake the decaying factor e^T^s(j-i)^ since A > 2. Taking also into account the factor due 
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j >' 

FIG. 4: Left: The ratio / {m?) between the typical momenta produced around zero and the average mass in every oscillation for 
the W (dashed blue line) and Z bosons (continuos red line) as a function of the number of oscillations. This ratio is significantly 
smaller than 1 for all crossings, which allows us to consider the produced gauge bosons as non-relativisitic. Right: Succesive 

spectral distributions k'^nk{l^)e'^^^''=^ at different j's, including the volume factor k^. One can see the predicted p29|) 
slow displacement of the maxima of the distribution. The x-axis is given in terms of a; = fc/(fc, (1)) 



to the expansion of the Universe, the first result we can read from here is that only for those values of j for which 
(j — l)logA — 21ogj > jF^{j — 1), the resonant effect will dominate over both the perturbative decay and the 
expansion rate. 

Note that inside the integral ()127p . the function u2g-(-D+s x;i=2 » ^^)" has a maximum at a value Up = {D + 



BEL. 



-2/3\-l/2 



which implies that that the typical (comoving) excited momentum is 



(129) 



In the right hand side of Figure 21 one can easily observe this behaviour: the value of the momentum at which the 



distribution peaks, slightly moves to smaller values, according to kp (x {D + B J2i=2 



,:-2/3)-l/2 



Thus, it is noticeable 



that the typical momentum k of the resonant fluctuations is always of order fc*(l), independently of how many 
oscillations the Higgs performs. The reason of this is that the parametric resonance effect builds up initially from the 
spectral distribution nfc(l+), which only filters k < fc*(l), see Eq. (jll7p and Figured! 

On the other hand, the ratio kp/{m)j for both W and Z bosons, between the typical momenta produced around 
zero and the average masses in every oscillation is shown in Fig. 2) In particular, it is easy to estimate the evolution 
in time of such a ratio, in terms of the resonant parameters (j85p . as 



{kp/uj 



« 1 Vj , 



(130) 



Taking into account that the previous ratio (|130p is a decreasing function with j , as well as its dependence with the 
gauge couplings g~^^^, we can conclude that the vector bosons produced at the bottom of the potential are always 



non-relativistic. This justifies a posteriori the calculation of the energy of the fermions as Epij) 



'jnz,w)j 



see 



Eq. ([M)) in section IIV CI Extrapolating Eq. (|130p to the case of fermions we realize that the produced particles at 
each zero crossing would be mainly relativistic, due to the smallness of the Yukawa couplings, being the only exception 
the production of top quarks. 

The energy density transferred to the fermions between the j-th and the [j + l)-th scatterings, will be 



^Pf[]) = 6 [(1 - e-'-^W)nz(j + )i?^^,(j) + 2(1 - e-''^^'^^^^)nw{t)EF„{]) 



((1 - e-''^^(^'))e-'^^^^(^-i) + 2cos 6*^(1 - e-'^«-^(^))e-'^"''^^(^-i)) , 



(131) 



where we have used the energy of the fermions (j99p and defined a momentum scale independent of the gauge couplings. 
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common to both bosonic species, as 



(132) 



The gauge coupHng dependence is indeed incorporated on the definition of the parameter 

(133) 



(cos 6*1^)3^5/2 (aKXond)^ 

which modulates again the strengh of the effect as e oc (/f. 
The total energy density transferred into the fermions will be 



8/3 

(134) 



8/3 



1=1 '' 
and the ratio of such an energy to that of the inflaton, 

M^j-^-^^^^^E^Mofi) . (135) 

with App{i) given by Eq. (|13ip . Here we can clearly see the two competing effects; that of the perturbative decay of 
the bosons, given by the factors of the form (1 — e^''-^'^^'>)e^^^^^^~^\ which tend to decrease the rate of production 
of bosons and fermions, while the factors e^'^^'' encoded in the form of the gaussian approximation, describe the 
resonant effect due to the accumulation of previously produced bosons and fermions. Initially, the perturbative decay 
will prevent the resonance to be effective. However, after a certain number of oscillations (a number that we will 
estimate next), the resonant effect will overtake the perturbative decays and parametric resonance will be developed 
as usual, as if the produced bosons would not decay perturbatively during each semi-oscillation. 

In order to estimate the time in which the perturbative decays stop blocking the parametric resonance effect, we 
can evaluate numerically when the expression e"^'"^^^"^' becomes subdominant versus e^'^ A^fcC*) , particular, 
we can evaluate the ratios, for either W 01 Z 

7rs(j - 1) 

for the fastest growing mode k = kp ()129|) . and find the number of semi-oscillations jn for which the previous ratio 
becomes greater than one, cr ^ 1. We find jj^ w 62 for the W bosons and jjf w 360 for the Z bosons. The fact that 
parametric resonance becomes important much earlier for W^s than for Z's is not a surprise, since their decay rate (j53|) 
differ in a factor ^z/jw ~ 2.4, which simply means that there are many more W bosons surviving per semioscillation 
than Z bosons. Therefore, the combined preheating of the W bosons is much faster driven into the parametric-like 
behaviour, while the evolution of the Z bosons is much more affected by the perturbative decays, delaying (or even 
completely preventing) the development of parametric resonance. Obviously, after a dozen of oscillations, the transfer 
of energy from the inflaton to the gauge bosons will be completely dominated by the channel into the W bosons, since 
by that moment they will be fully resonant while the Z bosons will still severely affected by their perturbative decay. 

Finally, to conclude this section and achieve an overall complete picture of all the details, let us also estimate the 
transfer of energy from the inflaton to the gauge bosons. In particular, the total energy transferred to them just after 
the j-th scattering, psij), is given simply by 



PB 



(j) = 3 {nzi3+){mz)j + 2nwU^){mw)j) , (137) 

where we have used the fact that the gauge bosons are non-relativistic and have 3 polarizations. Therefore, the ratio 
of the energy of the gauge bosons to the energy of the inflaton, can be expressed as 

esU) ^'-^-{j + iy i-^X f-~^^^^^^)^'-'^ (e--0) + 2cos.3,e---0)) , (138) 



where we have used Eq. (|127p and e, defined in Eq. (jl33p . modulates again the amplitude of this growing function. 
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A 


0.2 0.4 0.6 0.8 1.0 


Jeff 

AP) 

Jeff 


74 64 60 57 55 
79 69 64 61 59 



TABLE I: Number of semi-oscillations of the Higgs required, as a function of A, for an efficient transfer of energy from the 
inflaton to the gauge fields and/or to the fermions. 



Using Eqs. (|135p and p38p , we can estimate the time in which finaUy the energy of the inflaton would be transferred 
efficienty to the fermions or the bosons. Defining that moment, respectively, like £F(icft) = 1 and ^^(joff) = 1, 
one obtains the numbers in Table I. Note that the bosons receive the transfer of energy from the inflaton before 
the fermions, since by the time that parametric resonance overtakes the perturbative decay, the fraction of bosons 
decaying (per semi-oscilation) into fermions is very small and, therefore, the fraction of newly added fermions is less 
and less important, while the amount of produced bosons is more and more prominent. Note also that the number of 
oscillations jcff required for an efficient transfer of energy, depends on the parameter A, although the overall order of 
magnitude does not change appreciably. 

Unfortunately, as we will see in the next subsection, before reaching the stage in which ep.B ~ 1, the backreaction 
of the produced gauge fields into the homogeneous Higgs condensate will become significant, and it will have to be 
taken into account. 



VI. BACKREACTION 



Let's now calculate the backreaction from the W and Z bosons into the Higgs condensate. Neglecting the vectorial 
nature of the bosons, the effective equation of the Higgs condensate can be written, in the Hartree approximation, as 



4C Xdx 



x = o. 



(139) 



where there is a field for each polarization of each gauge boson species, such that ~ g\ and — / cos^ Ow, as 
usual, for the W and Z bosons, respectively. From here, performing the derivative, one obtains for the effective Higgs 
frequency 



u;2 ^ + 



ag'^Mj 



(140) 



where the second term in the r.h.s. should be summed over polarizations and species. For the fraction of time of each 
semi-oscillation, during which the Higgs frequency evolves adiabatically, we can use the correlation function 



with (pk{t) expressed as 



y2ZJ(fc) vMfc) 

Thus, one can compute the expectation value of the bosonic fields (components) 



(141) 



(142) 



27r2a3 



dkk'^\ipk\ 



27r2a3 



dkk^ n 



1 1 
27r2a3 gMp ^\ _ e-"«lxl 



1 -I- cos 



Arg ak + Arg 



(143) 



where, to obtain the last exppresion we have used ujk = y^-Vl — e^"'*lxl , = IctfcP — 1 = rife and J u!k{t')dt' 



(7r/M)X;"=i(w>j, with = K dt'Lo{t'). Following [Si], since we don't know the accummulated phases of afe 



24 



A 


0.2 0.4 0.6 0.8 1.0 


J/backr 


67 57 52 50 48 



TABLE II: Number of semi-oscillations of the Higgs required, as a function of A, for the backreaction of the gauge fields into 
the Higgs background to become significant. 



and Pk, we will write 



{>n^) ^ ^ U^A en. Y {'■'),)] . (144) 

' gMp Vl - e-"«lxl L \M "/J' ^ ' 



with A <1 and = (27r^a^) ^ / dkk'^nk- 

From here, one can define the effective frequency of the Higgs condensate as 



(145) 



The backreaction of the gauge boson fields over the Higgs field, will be non-negligible when the last term in the r.h.s. 
of the previous expression becomes of the order of M^. In terms of the number densities of the Z and W bosons, i.e. 
summing the contribution over polarizations and species of all the fields that back react, this will happen at a time 

= jn/M, 

Backreaction ^ (n,{j)/ cosO^ + 2n^U)) > ^^IxfeOK^^IxfeOD^^^M^ ^ 
V / oag2 



where we have expanded V e^^'^l^l — e°'^\^\ ~ (aK|x|)^^^, which is certainly accurated after a couple of dozens of 
ocillation, since |x(ira)l oc i. Substituying the averaged value per semi-oscillation x{t) (x(^))jj then we will take 
aK|x(ij)| ""^^J"'' s'm{x)) = ^ ""^^J"'' . Using the analytical expressions for the occupation numbers (|127p we 

can translate the above condition into the following one 

Thus, if we find numerically the number of semi-oscillations of the Higgs, jbackr7 for which j > jbackr fuUfiUs the 
above condition, then we know the moment in which backreaction of the bosonics fields becomes significant, tbackr ~ 
TTjbackr/M. Notc that the above condition depends on A both in the left- and right-hand sides. In particular, the A 
dependence is rather weak in the constant of the right-hand side of the inequality, since it goes as A^/*, while in the 
left-hand side, it enters through the exponentials so it can change the number ribackr in a more significant manner. 
Taking values of A between 0.2 and 1.0, we obtain the numbers in Table II. 

We clearly see that backreaction seems to become important at a time slightly earlier than that at which we 
were expecting the Higgs to have transferred efficiently its energy to the bosons and fermions. This means that our 
analytical estimates of these transfers were biased, and a careful numerical study of the process is required. Beyond 
backreaction, the strength of the resonance very quickly decreases due to the increased frequency of oscillations of the 
Higgs. Eventually, the broad resonance driving the production of gauge bosons and thus their decay into SM particles 
becomes a narrow resonance and finally shuts off. From then on, the inflaton will oscillate like a matter field while 
the produced particle will redshift as radiation, its effect on the expansion becoming negligible after a few hundred 
oscillations. 



VII. CONCLUSIONS 



We have studied the different stages of reheating after inflation in a model where the role of the infiaton is played 
by the Higgs field of the Standard Model of particle physics with a non-minimal coupling to gravity. Infiation in 
this model takes place at the GUT scale, along the lines of the Starobinsky model of infiation since a conformal 
transformation makes these two models indistinguishable from the point of view of infiation. The usual difficulty with 
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large self couplings of the Higgs is tamed here by the inclusion of a large non- minimal coupling to gravity, ^ ^ 10^, 
which nevertheless does not leave any signature at low (electroweak) scales due to the fact that the Higgs field acquires 
a vacuum expectation value and does not evolve at present, while the local spacetime curvature is negligible. 

The advantage of this model of inflation for the study or reheating after inflation is that all the couplings of 
the Higgs-Inflaton to matter fields are known at the electroweak scale, and can be extrapolated to the GUT scale 
using the renormalization group equations, and therefore one can study in detail the process of reheating of the 
Universe, without having to impose ad hoc assumptions about their values. The surprise is that the process becomes 
more complicated than expected, and a series of subsequent stages take place, where essentially all different types of 
particle production mechanisms at preheating occur. Moreover, since the Standard Model couplings of the Higgs to 
gauge and matter fields are non-negligible, nor are their couplings among themselves, the process of non-perturbative 
decay via parametric resonance is mixed with the usual perturbative decays of the decay products, which complicates 
things significantly. 

Inflation ends at values of the Higgs field of order the Planck scale and goes through a brief stage of tachyonic 
preheating soon after the end of inflation. The passage is so short that particle production is not significant at that 
stage. The same occurs with the production at the inflection point. Finally the Higgs-Inflaton field starts oscillating 
around the minimum of its potential with a curvature scale of order 10^^ GeV. At this stage, particle production 
occurs whenever the Higgs passes through zero, creating mostly vector gauge bosons W and Z. These gauge bosons 
acquire a large mass while the Higgs increases towards maximum amplitude and start to decay into all Standard 
Model leptons and quarks within half a Higgs oscillation, rapidly depleting the occupation numbers of gauge bosons, 
like in instant preheating. However, the fraction of energy of the Higgs that goes into SM particles is still very small 
compared with the energy in the oscillations, and therefore the non-perturbative decay is slow. This implies that 
a relatively large number of oscillations take place before a significant amount of energy is transferred to the gauge 
bosons and fermions. 

The amplitude of Higgs boson oscillations decreases as the Universe expands in a matter-like dominated stage with 
zero pressure. Eventually, this amplitude is small enough that the gauge boson masses are not large enough for 
inducing a quick decay of the gauge bosons and these start to build up their occupation numbers very rapidly via 
parametric amplification. The question whether this effect can give rise to the production of a significant Gravitational 
Wave Background (GWB) potentially observable today remains to be addressed. Several papers have studied recently 
such an issue in the chaotic and hybrid models of inflation [4l| . but in the present model, we don't have simply a 
parametric resonance phenomena but a combined preheating effect which, perhaps, could modify the properties of such 
a GWB. Similar arguments would affect also the production of magnetic fields at preheating [42| or even electroweak 
baryogenesis [isj . 

After about a hundred oscillations the gauge bosons produced backreact on the Higgs field and the resonant 
production of particles stops. The Higgs field acquires a large mass via its interaction with the gauge condensate and 
preheating ends. From there on, both Higgses and gauge fields decay perturbatively until their energy is transferred 
to SM particles. Since the stage after backreaction is very non-linear and non-perturbative, it cannot be solved 
analytically and we have to resort to numerical studies in the lattice. We leave the description of our numerical 
studies to a future publication. 

Note added: Upon completion of this paper, we received through the arXiv the preprint of Bezrukov et al. [i^ j. 
where they also study preheating in the uMSM. Although the formalism is common to both, our conclusions are 



somewhat different from those of Ref. [4j|. We find that the Higgs decay into gauge bosons is significantly faster, 
and that backreaction occurs much before thermalization. We thus think it is not possible to determine the reheating 
temperature without a careful numerical analysis with lattice simulations. 

A few days after this work was completed, Ref. [45] was also posted in the arXiv, performing an analysis of the 
2-loop quantum corrections to the running of all the parameters involved in the model. This paper allows for a 
different range of the Higgs self-coupling, which is compatible with the range (j38p although more restricted. The main 
result of Ref. ^45] is a relationship between the Higgs mass and the spectral index which, in princi ple, could be tested 
in the future against data from PLANCK and the LHC. Similar conclusions where found in Ref. [iq . 
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